Gold Clusters Sliding on Graphite: 
a Possible Quartz Crystal Microbalance Experiment? 



o 
o 

(N 



S. Pisov, 1 - 2 E. Tosatti, 3 - 1 U. Tartaglino, 3 and A. Vanossi 4 

1 International Center for Theoretical Physics (ICTP), P.O.Box 586, 1-34014 Trieste, Italy 
2 University of Sofia, Department of Atomic Physics, 5 James Bourchier Blvd., 1164 Sofia, Bulgaria 
3 International School for Advanced Studies (SISSA) and Democritos, Via Beirut 2-4, 1-34014 Trieste, Italy 
J * CNR-INFM National Research Center S3 and Department of Physics, 
University of Modena and Reggio Emilia, Via Campi 213/ A, 41100 Modena, Italy 

(Dated: December 4, 2006) 

A large measured 2D diffusion coefficient of gold nanoclusters on graphite has been known ex- 
perimentally and theoretically for about a decade. When subjected to a lateral force, these clusters 
should slide with an amount of friction that could be measured. We examine the hypothetical 
possibility to measure by Quartz Crystal Microbalance (QCM) the phononic sliding friction of gold 
clusters in the size range around 250 atoms on a graphite substrate between 300 and 600 K. As- 
suming the validity of Einstein's relations of ordinary Brownian motion and making use of the 
experimentally available activated behavior of the diffusion coefficients, we can predict the sliding 
friction and slip times as a function of temperature. It is found that a prototypical deposited gold 
cluster could yield slip times in the standard measurable size of 10 _9 sec for temperatures around 
450 — 500 K, or 200 C. Since gold nanoclusters may also melt around these temperatures, QCM 
would offer the additional chance to observe this phenomenon through a frictional change. 
Prepr int of: J. Phys.: Condens. Matter 19 (2007) 305015. 
Link: |http://stacks.iop.org/0953-8984/19/305015| 
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I. INTRODUCTION 

Understanding the diffusion mechanisms and the fric- 
tional properties of aggregates of atoms or molecules of 
nanometric size (nanoclusters) on surfaces is important 
from both fundamental and technological view points. 
Growth of new materials with tailored features, as, e.g., a 
structure controlled down to the nanometer scale, is one 
of the active research fields in physical science. Differ- 
ent experimental techniques can be used to build nanos- 
tructured systems; however the requisites of control (in 
terms of characterization and flexibility) and efficiency 
(in terms of quantity of matter obtained per second) 
are generally incompatible. The main advantage of the 
cluster-deposition technique is that one can carefully con- 
trol the building block (i.e., the cluster) and character- 
ize the formation processes. The behavior of these de- 
posited nano-objects is often distinctive (mostly due to 
their large surface-to- volume ratio), being qualitatively 
different from those of their constituent parts and from 
those of bulk material. In particular, they may present 
properties that vary dramatically with size. 

Under different circumstances, various individual 
mechanisms (single-atom-like processes) can be respon- 
sible for cluster motion, such as evaporation and/or con- 
densation, diffusion of particles along the cluster edge, 
and motion of misfit dislocations. Differently, one of the 
remarkable experimental observations of the last decade 
is concerned with gliding- like dynamics of compact solid 
gold clusters as a whole (concerted jumps with conserva- 
tion of size and shape). 

After room temperature deposition, gold clusters, com- 



prising typically 250 atoms (radius of ~ 10 A), have 
been repeatedly observed to diffuse on Highly Ordered 
Pyrolytic Graphite (HOPG) surfaces with surprisingly 
large, thermally activated diffusion coefficients in the 
range 10~ 7 — 10~ 5 cm 2 /sec already at room temperature 
[l[; a similar behavior was reported for Sb23oo clusters. 
A large mobility was recently observed also for activated 
surface diffusion of close-packed hexagonal clusters Irj 
and Irig 0. The detailed atomic mechanisms for diffu- 
sion of gold clusters, on which we focus here, was stud- 
ied theoretically and simulated by Molecular Dynamics 
(MD) by Luedtke et al. [3|, who pointed out the coex- 
istence of short and long jumps, the latter assimilating 
the process to a Levy flight. Later work by Lewis et 
al. [J] and by Maruyama [5j further explored the cluster 
diffusion mechanisms. 

The observable diffusion of sizeable clusters as a whole 
raises an interesting question in the context of nanofric- 
tion. If one could manage to apply a sufficient lateral 
driving force to the clusters, they would drift under its 
action. By the fluctuation-dissipation theorem, the lat- 
eral cluster drift mobility should then be related (leav- 
ing Levy flights aside for the time being Q) to diffusion 
through Einstein relation. The measurable diffusivity of 
deposited clusters could thus lead to a measurable fric- 
tional dissipation in experimental apparatuses such as a 
Quartz Crystal Microbalance (QCM) This is pre- 

cisely the issue we propose to explore synthetically in this 
paper. 

Specifically, we wish to examine a hypothetical, yet 
very practical, QCM experimental case where sliding 
friction could be measured for gold clusters adsorbed 
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on graphite, for sizes around 250 atoms and tempera- 
tures between room temperature and about 600 K. As 
it will turn out, the experimentally observed temper- 
ature dependence of the diffusion coefficients (obeying 
an Arrhenius-type activated law) predicts, through Ein- 
stein's relation, easily measurable QCM slip times at 
temperatures of the order of 450-500 K, possibly even 
lower, which anyway appear well within the reach of ex- 
isting QCM setups. 

An additional interesting effect is that small gold clus- 
ters will tend to "premelt" already at much lower temper- 
atures than the bulk T m = 1336 K. Buffat and Borel Q 
analyzed free gold clusters for premelting. Their data ac- 
tually suggest melting not too far from room temperature 
for sizes around AU250 or thereabout. In addition one 
may expect the solid to acquire some additional stabil- 
ity against the liquid when adsorbed on the flat graphite 
substrate. So the actual temperature where the adsorbed 
Au cluster melts into a (partly) wetting Au droplet is 
difficult to predict and will have to await experimental 
scrutiny. Our point however is precisely that a sudden 
change of QCM sliding friction as T increases should be 
generically observable, and will signal that the cluster 
has melted. 



II. MODELLING 

The simplest theoretical model for cluster diffusion 
could be a model of jump over simple monatomic steps. 
A single jump mechanism predicts diffusion coefficients 
given by 

where a is the jump step, v is the attempt frequency, 
and E^ is the energy barrier to be overcome. How- 
ever, even if v is taken as large as the Debye frequency 
(y rs 10 12 s _1 ) , with a ~ 0.3 nm, the prefactor is 10~ 3 
cm 2 /sec [Toj . many orders of magnitude smaller than ex- 
perimental value of 10 3 cm 2 /sec [l|. 

Luedtke et al. Q simulated gold cluster diffusion nu- 
merically using a many-body embedded atom (EAM) po- 
tential for the interaction among the gold atoms, and a 
two-body Lennard- Jones (LJ) potential for the interac- 
tion between gold and carbon atoms. The LJ poten- 
tial parameters were cau-c — 0.01273 eV and (tau-c = 
2.9943 A. The Au — C LJ potential was given different 
parameters cau-c — 0.022 eV and oau-c = 2.74 A in 
a later simulation by Lewis et al. Both works un- 
covered a cluster diffusion mechanism that has a very 
interesting nature. There are long "sticking" periods of 
short-range local dwelling without too much diffusion of 
the cluster. Every now and then, there appear rare but 
important events where the cluster "slips" - it actually 
glides - over relatively long distances. The sticking plau- 
sibly corresponds to trapped cluster configurations, occa- 



sionally abandoned thanks to large and rare fluctuations. 
During one such fluctuation the two crystalline surfaces 
in contact, the gold cluster and the graphite plane, be- 
have as hard incommensurate sliders [ll| . and thus slip 
for a while essentially freely. 

This reasoning is supported by existing numerical re- 
sults based on molecular dynamics simulations. Many- 
particle clusters that are incommensurate with the sub- 
strate were shown [12J to exhibit very rapid diffusion, 
with their paths akin to a Brownian motion induced by 
the internal vibrations of the clusters and/or the vibra- 
tions of the substrate. The behavior of an incommensu- 
rate object moving "as a whole" is in striking contrast 
with other diffusion mechanisms, especially observed for 
clusters epitaxially oriented on the surface, where the 
cluster motion results from a combination of single-atom 
processes (e.g., evaporation, condensation, edge diffu- 
sion, etc.). The latter mechanisms, giving rise to rel- 
atively low diffusion coefficients (D ~ 10~ 17 cm 2 /sec), 
are likely not significant in cases where the mismatch be- 
tween cluster and substrate is large and/or their mutual 
interaction is weak. In this view, a gold cluster adsorbed 
with a (111) face onto a graphite substrate forms an ideal 
mismatched (incommensurate) system with two very stiff 
mating surfaces that interact only weakly, thus represent- 
ing a desirable sliding situation. 

As shown by the simulation work by Lewis et al. 0], 
this stick-slip type dynamics still leads - despite its large 
difference from the simple jump model over monatomic 
barriers - to an activated diffusion coefficient of the form 
[T] above, although now with different coefficients from 
monatomic diffusion. That work however also clarifies 
the quantitative inadequacy of the LJ interaction model, 
which for a 250 atom cluster leads to a prefactor of about 
2 x 10~ 3 cm 2 /sec and an effective barrier of 0.17 eV, 
against experimental values of 10 3 cm 2 /sec and 0.5 eV 
[l| . The Arrhenius activated behavior fit to these experi- 
mental data is shown in Fig. [I] The discrepancy between 
simulation and experiment was highlighted by Maruyama 

The source of discrepancy between simulations and ex- 
periment is probably related to a large error in the ad- 
sorption energy. An experimental estimate for adsorption 
energy of Au on graphite is E a = 0.64 eV [131 ] . ab-initio 
calculations suggested even l arg er values for the adsorp- 
tion energy of about 0.9 eV [14|. On the other hand the 
LJ model potentials imply a much smaller adsorption en- 
ergies E a ~ 0.15 eV. With this smaller adsorbtion energy, 
a cluster could experience a smaller effective energy bar- 
rier when trying to disentangle itself from the graphite 
substrate. Thus the adsorption energy error, at first sight 
irrelevant, might we suspect be at the origin of the large 
discrepancy of effective barriers in cluster diffusion (0.5 
eV against 0.17 eV in Lewis's LJ modeling, and probably 
even smaller in Luedtke's). 

While it will be important to fix these problems in 
future simulations, we can for the time being content 
ourselves with this qualitative understanding, and simply 
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tion coefficient in ([3]), we obtain for the slip time 

r = Dm Au N/k B T. (4) 

Through this relation, assuming approximately N = 
50 atoms directly touching the graphite substrate and, 
from the experimental data of Bardotti et al. lj, a diffu- 
sion coefficient D ~ 10~ 5 cm 2 /sec for AU250 nanoclusters 
at 300 K, we get a rough estimate of the room temper- 
ature slip time r ~ 10~ 12 sec. While this value is at 
least two orders of magnitude too small to be observed 
by QCM, the situation, as discussed below, can change 
drastically upon heating. 



IV. A POSSIBLE QCM EXPERIMENT 



FIG. 1: Temperature dependence of the diffusion coefficient 
of AU250 on graphite. The values of Do and Ed are chosen 
to fit the experimental data by Bardotti et al. [1J with an 
Arrhenius-type law. 

discuss cluster diffusion as a regular Brownian diffusion, 
although with a rather unusual activation mechanism. 

III. ESTIMATE OF THE SLIP TIME THROUGH 
EINSTEIN'S RELATION 

When an adsorbate island or a cluster is forced to slide 
on a substrate, the dissipation of energy due to kinetic 
friction can be characterized by a specific tribological 
quantity, the slip time r, defined as the time taken by 
the cluster speed, initially set to be nonzero, to drop to 
1/e of its original value. This is a relaxation time as- 
sociated with the cluster-momentum fluctuation, and is 
connected to the interfacial friction coefficient i] through 
the relation r) = p/r, where p is the mass per unit area 
of the cluster. Thus, defining p = rriAuN/A (N denoting 
here the number of gold atoms of the cluster in direct 
contact with the graphite substrate over an area A) , the 
interfacial friction coefficient can be rewritten as 

77 = m Au N/AT. (2) 

Our assumption that cluster diffusion enjoys the prop- 
erties of Brownian motion implies the applicability of the 
fluctuation-dissipation theorem, which in turn causes the 
force free diffusion and the friction under a sliding force 
to be connected through Einstein's relation 

DrjA = k B T. (3) 
By substituting expression Q for the interfacial fric- 



The problem of how to measure interfacial friction in a 
quantitative manner has remained long unsolved, while 
knowledge of this property would provide important in- 
formation relevant to a wide variety of problems. Sliding 
friction of adsorbed films are very effectively probed by 
QCM. In a QCM experiment, a substrate is laterally os- 
cillated with typical frequencies / ~ f 5 MHz and ampli- 
tudes of a ~ 100 A. This exerts on a iV-atom adsorbate 
island to a typical inertial force F = N(2tt) 2 M f 2 a, which 
when sufficiently large can depin the island and make it 
slide. The sliding frictional dissipation lowers the oscil- 
lator's quality factor Q, and the result is conventionally 
measured by the slip time r = d(Q~ 1 )/df. As was said 
earlier, the interfacial friction coefficient may then be de- 
rived from Eq. [2j Sliding friction vanishes in both limits 
t = (film locked to the substrate), and r = oo (super- 
fluid film), and typical measurable slip times are in the 
order of 10" 8 - f0" 9 sec 

For low speed, hard sliders on hard flat substrates far 
from their melting points, there is no wear, and all fric- 
tional dissipation occurs via phonon creation (in met- 
als, also creation of electron hole pairs). We concen- 
trate here on phononic dissipation, whose origin stems 
from the oscillatory potential energy variation associated 
with lateral cluster displacement, the so-called "corruga- 
tion" . A recent analysis of solid xenon monolayers sliding 
over different substrates suggested for example that the 
viscous friction coefficient 77 increases quadratically with 
corrugation [l6| 

V = Vsubs + bUg (5) 

where i] su b s is the dissipation arising from both phononic 
and electronic friction, b is a coefficient which depends 
weakly from the substrate surface, and Uo is the potential 
corrugation amplitude between adsorbate layer and sub- 
strate, controlling the phononic dissipation. Experimen- 
tal upper bounds for r) su bs &re 0.08 ns" 1 (forXe/Cu(m), 
Xe/Ni(lll), and Xe/graphite). 
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FIG. 2: Predicted temperature dependence of the slip time of 
a AU250 on graphite, r is evaluated considering an interfacial 
contact area of TV = 50 Au atoms. 

We propose here a possible QCM experiment to study 
the frictional dissipation of gold nanoclusters sliding on 
a graphite substrate heated above room temperature. 
Raising temperature will bring about QCM-measurablc 
slip times. Working at temperature even as high as ~ 600 
K does not affect the structural properties of graphite, 
and should be entirely possible. Assuming for the mo- 
ment the Au clusters to remain solid, we can simply ex- 
trapolate at higher temperatures the diffusion and slip 
times as done in Fig. O The predicted cluster slip time 
should for a 250 atom cluster reach QCM measurable val- 
ues of 10 -9 sec between T = 450 and 500 K, where the 
corresponding diffusion coefficient D ~ 10 -2 cm 2 /sec. 

V. AU CLUSTERS MELTING AND 
NANOFRICTION 

As was recalled in the introduction, small clusters of 
radius R have a premelting temperature T m (R) well be- 
low the bulk melting temperature of the same material. 
Figure [3] reproduces the original data by Buffat and Borel 
@ , showing that a small gold cluster of diameter ~ 20 A 
will generally melt at room temperature or below. Pre- 
melting of gold clusters has also been extensively simu- 
lated, notably by Ercolessi et aLp"^ 

The premelting temperature data for free gold clusters 
are well fit by the Gibbsian formula 

T m (R)^T m (oo)(l--^A 7 j, (6) 

with T m (oo) = 1336 K, and A7 = j sv - (j SL + ~f LV ) re 



600 mJ/m 2 . Here, pi is the density of the liquid, L is 
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FIG. 3: Experimental values of the melting point temperature 
of free gold clusters as a function of their diameter. From 
Ref. 0. 



the latent heat of melting per unit mass, and the 7's 
denote the free energies of the three solid- vapor (SV), 
solid- liquid (SL), and liquid- vapor (LV) interfaces. 

Even though a deposited cluster might generally be 
more stable than a free one due to the stabilizing effect 
of a hard substrate, it will still premelt very readily upon 
heating. At melting, one should expect the nanoclus- 
ter friction to change very significantly. In fact the nan- 
oclustcr will turn into a droplet, with some characteristic 
wetting angle, and obviously losing its original rigidity. 
The droplet dissipation should be dominated by viscosity, 
with a damping whose magnitude is at this stage difficult 
to anticipate. The jump between rigid cluster sliding and 
viscous droplet dissipation should nonetheless be a strong 
feature, readily observable by QCM, taking place below 
or near 600 K, depending on size, for clusters below 20 
Ain diameter. 

In summary, we predict that QCM measurement of 
small gold clusters on a graphite substrate should lead 
to observable inertial sliding friction with measurable 
slip times at temperatures of a couple hundred centi- 
grades. The relatively sudden melting of deposited clus- 
ters should also become measurable in the form of a sud- 
den jump of slip time upon heating, at temperatures that 
should depend very strongly on cluster size. 
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